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The limiting bubble wall velocity during a first-order electroweak phase transition is 
of interest in scenarios for electroweak baryogenesis. Khlebnikov has recently proposed 
an interesting method for computing this velocity based on the fluctuation-dissipation 
theorem. I demonstrate that at one-loop order this method is identical to simple, earlier 
techniques for computing the wall velocity based on computing the friction from particles 
reflecting off or transmitting through the wall in the ideal gas limit. 
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There has recently been considerable interest in the limiting bubble wall velocity 
during a first-order electroweak phase transition, era The dynamics of the bubble wall 
bears on scenarios for electroweak baryogenesis mediated by the anomalous violation of 
baryon number in the Standard Model. §h0 During the phase transition, bubbles of 
the symmetry-broken phase are nucleated inside the symmetric phase, and these bubbles 
then expand to convert the entire Universe into the symmetry-broken phase. At zero 
temperature, in vacuum, the velocity of the bubble walls would asymptotically approach 
the speed of light; at finite temperature, however, collisions of particles off of the bubble 
wall provide a source of friction and force a smaller limiting velocity. A simple model 
for this friction is to compute the momentum transfer to the wall from an ideal gas of 
thermally-distributed particles colliding with, and either reflecting from or transmitting 
through, the bubble wall.* This model is valid provided the thermal mean free path 
of those particles is large compared to the thickness of the bubble wall. (This is not a 
particularly good assumption for the real case of the electroweak phase transition, but 
is often used to establish a lower limit on the bubble wall velocity. ) Khlebnikovi has 
recently proposed an interesting alternative method for computing the bubble wall velocity 
u, based on the fluctuation-dissipation theorem, for cases where u C 1. When he puts 
his general formalism into practice with a sample one-loop calculation, he claims to find 
slightly different results than expected from other methods. In this paper, however, I shall 
show that the one-loop result obtained from the fluctuation-dissipation theorem is identical 
to the simple model of reflection and transmission in the limit of small wall velocities and 
large mean free path. (Khlebnikov's general formalism is still of interest because it applies 
to higher- loop calculations as well, always provided « <C 1. For example, the effects of 
rescattering, which go beyond the approximations considered in this paper, originate in 
higher-order loops.) 



This has been discussed, in various limits, in refs. flljp, 11 ,12|. 



I shall also reexamine a toy model calculation considered by Khlebnikov, where he 
found a logarithmic suppression of the bubble wall velocity. It will be possible to see the 
physical origin of this logarithm (and of another logarithm that he accidentally omitted) 
by using the picture of reflection and transmission at the bubble wall. 

1. Review of simple reflection/transmission model 

Figure [l] shows the qualitative form of the effective, finite-temperature Higgs potential 
when the phase transition begins: the symmetric = phase (phase 1) is unstable to decay 
into the asymmetric vacuum (phase 2). Figure |] shows qualitatively the profile of a bubble 
wall interpolating between the phases and which, as drawn, will move from right to left. 
Since particle masses depend on the value of (p, masses will be different in the two phases. 

If the mean free path is large compared to the wall thickness, then we may separately 
treat the collision of each particle in the plasma with the bubble wall. In the wall rest 
frame, the wall will be struck by a Lorentz-boosted thermal distribution of particles. I 
shall focus on the case of bosons. The usual Bose distribution n p = l/[exp((3E p ) — 1] is 
the form, in the thermal rest frame, of the Lorentz scalar l/[exp(/3p • v) — 1] where v is the 
four-velocity of the thermal bath. In the rest frame of the bubble wall, moving with wall 
velocity u relative to the thermal bath, this distribution is 



Take u to be in the z direction, and consider a particle of energy E and momentum 
p = (p±,p z ) incident from the symmetric phase. Both E and p± will be conserved. The 
momentum transferred to the wall will be 2pi z if the particle is reflected and p± z — p2 Z if 
it is transmitted, where 
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(1.1) 
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and mi and rri2 are the masses of the particle in the two phases. The pressure on the wall 
from particles incident from the symmetric phase is then 



/d pi 
j^^v lz 6(v lz )f pl [2pi z 7£(pi) + (pu - P2 Z )T(p 1 )] , (1.3) 

where lZ(pi) and T(pi) are the reflection and transmission coefficients for particles en- 
countering the wall and v% is the incident particle velocity pi/E. The pressure P2 from 
particles incident from the asymmetric phase is similar. t Also, the difference of the Higgs 
potential between the two phases exerts pressure on the wall. The net pressure is 

P net =P 1 -P 2 - AV ch (1.4) 

where V c \ is the classical, zero-temperature Higgs potential. (The finite-temperature effec- 
tive potential will appear momentarily.) 

Consider this result in the limit of zero wall velocity u, and change integration variables 
in ( |1.3|) from p z to E: 



P„ t(M= o) = / d '£)f ^ t {iW>(-A) (1 5) 

+ 2(p lz -p 2z ){Tl(p ± ,E)+T(p ± ,E)]0(plJ}-AV ch 

where I have used the fact from one-dimensional quantum mechanics that reflection and 
transmission coefficients do not depend from which side the particle approaches. 9(±P2 Z ) 
above is a short-hand notation for 9[±(E 2 —p\ — m^)], which specifies whether or not the 
particle is classically allowed in phase 2. The first term in braces above corresponds to 
particles from phase 1 which are not energetic enough to enter phase 2 and so must reflect 
from the bubble wall; the second term corresponds to more energetic particles, which may 



' I shall assume throughout that there the temperature does not vary appreciably across the 



wall. This is the case for the electroweak phase transition; see refs. Q and [12]. 



come from phase 1 or phase 2 and may reflect or transmit. Now use 1Z + T = 1 and 
integrate by parts to get 



P net (u = 0) = Tj ^ [In (1 - e-M) - In (1 - e'^)} - AV ch 



;i.6) 



which is simply the difference —AV e g in the one-loop finite-temperature effective potential 
between the two phases. @ 

At non-zero wall velocity, the net pressure ( |1 . 3| ) may then be written as 

d 3 pi 



Pnet = j ^^V lz 9(v lz )Af pl [2pi z K(pi) + (p lz -P2b)T(Pi)] 

- J ^^v 2z 6(-v 2z )Af p2 [2p 2z Tl(p 1 ) + (p 2z -Pi a )T(pi)] - A14 ff , 



;i-7) 



where Af p = f p — n p . In the limit of small wall velocity u, f p — n p ~ u(3p z n p (n p + 1), and 
( pT7| ) simplifies to 

P net = -AV eS +(3u J ^^ n (E)[n(E) + l]{2&e{-&) 



in 2 ) 



+ [2(p? a + pL)K(p ± , E) + (p lz - p 2z ) 2 T(p ± ,E)} 9(p 2 2z ) } + 0(i 

(1.8) 



where p\ z and p 2z will henceforth always be understood to be positive, according to ( |1.2|) . 
The limiting velocity of the wall is determined by setting P net = 0. 

Equation ( |1.8| ) is the formula that I shall show is equivalent to the one-loop fluctuation- 
dissipation result of Khlebnikov. S For the sake of simplicity and concreteness, I shall work 
with a generalization of the toy scalar model considered by Khlebnikov. The Lagrangian 
is 

C = ±(d<P) 2 - V{4>) + \m 2 - \M 2 m\ (1-9) 

where <p is the usual Higgs field, and $ is another particle in the theory that, analogous 
to W and Z bosons, has a mass that depends on the VEV of <p. [Khlebnikov's specific 

4 



model was M 2 ((f)) = Mq + 2gcf>.] In the presence of a background field <f> c i(z) describing 
the bubble wall in its rest frame, the $ field has the effective Lagrangian 

£ eS = ±(d<5>) 2 -±M 2 (Mz))$ 2 - (1-10) 

This is just the problem of a particle propagating in a classical potential. The S matrix 
may be found by solving the classical equation of motion, taking a solution with a single 
incoming wave, and extracting the amplitudes of outgoing waves. Fourier transforming 
(t, x±) to (E,p±) and noting that the latter are conserved, the equation of motion is just 

[-d 2 + M 2 (4> c \(z))\ d> = (E 2 - pi)$, (1.11) 

which is nothing more than the Schrodinger equation for a unit mass particle in the po- 
tential 2M 2 (^ c i(z)) with energy 2{E 2 — PjJ- For incident waves from the right or left, it 
has solutions shown qualitatively in fig. |||, and the reflection and transmission coefficients 



extracted from these solutions are the 1Z and T which appear in (|1.8|). 



2. Fluctuation-dissipation theorem at one loop. 

I shall now show that one also obtains the result ( |1.8| ) from the fluctuation-dissipation 
theorem at one-loop. Khlebnikov treats the difference in vacuum energies in fig. |asa small 
perturbation. Imagine turning off that perturbation, so that the vacua are degenerate. A 
bubble wall in this situation would then remain stationary. Now turn the perturbation 
back on. The non-equilibrium response of the system to this perturbation can then be 
related, via the fluctuation-dissipation theorem, to the equilibrium calculation of a retarded 
correlator. Specifically, assuming the bubble wall velocity is small, he finds 

u = AV e ff lim 
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dzdz' d z <fi c i(z) Imn R (w; fcj_ = 0; z, z') d z >(j) c i(z') 



(2.1) 



where ITr is the retarded self-energy of the Higgs field <f>. ITr is projected above onto the 
translational zero-mode d z (j> c \{z) of the bubble wall. In position space, ITr = ([(f)(r), 4>{r')]) 
is a function of two space-time coordinates. I have Fourier transformed the t, x, and y 
components above to to and k± but have left the z coordinates in position space. The trick 
is to now compute the desired limit of ITr in the background of the bubble wall (f> c \(z). I 
shall examine this problem at one-loop order. 

For simplicity, return to the scalar model ( |1.9|) and focus on the one-loop contributions 
to ITr from $, which are shown in fig. |j. Only fig. |]b has an imaginary part. To get 
its contribution to IulIEr, first consider evaluating the self-energy in Euclidean space for 
Euclidean frequencies iv. Figure [|b gives 




( 



(p + iv) 2 -p\ 



1 



A 2 



) 



(2.2) 



where the sum is over discrete po = i2imT, A 2 is the operator 



A 2 



d 2 z +M 2 (<p cl (z)), 



(2.3) 



and 



M'(z) =d 4> M\(t )cl {z)). 



(2.4) 



Now let fj^ (z) be the eigenfunctions of A 2 defined by 





(2.5) 



and normalized to 




(2.6) 
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Later on, we shall see that a useful basis for the f ± will be the solutions depicted in fig. [| 
representing scattering of particles incident from a definite side of the wall. By inserting a 
complete set of states, ( |2.2|) may be rewritten as 



TL{iu;0;z,z') = \m \z)M \z') [ dn 2 dn' 2 Tj2 / 

2 J Po J 



d 2 p 



_L 



1 



* f n \Z) o "> > • 

Po-P±- K 



27T 



(p + iv) 2 — p\ — K 



f b Az), 

(2.7) 



where sums over a, b = ± are implicit. At this point, the steps to obtain ImilR are exactly 
the same as they would be in the case of free scalar propagators and may be found, for 
example, in refs. |13|J^|. First the frequency sum over po is performed by the usual contour 
trick to obtain 



U(iu;0;z,z') = -M'(z)M'(z') / d«W7r (z')fZ(z') 



n(E) 

2E 
n{E') 

2E' 



1 



+ 



1 



(E + iv) 2 - (p 2 ± + k' 2 ) (E - iv) 2 - (p 2 ± + k' 2 ) 

1 1 

+ 



+ (T-independent terms), 



where 



p^ + K 2 , E' = Jp\+K' 2 



(2- 



(2.9) 



The retarded self-energy is then obtained by replacing iv by uj—ie. It is now straightforward 
to take the imaginary part, 

d 2 p± 



lmU R (u;0;z,z') = lM'(z)M'(z') j dK 2 dK' 2 fT V)tfCO j 



(27TY 



TV 



AEE 



7 [n{E) - n(E')} [5(u + E — E') — 8{uj — E + E') 



(2.10) 



and to take the 10 — > limit. Switching integration variables from k to E, 



lmU R (u;0; Zl z') ^ rr^M'(z)M'(z') J |^ J 



dE 



[2.11) 



x n(E)[n(E) + (*)/*(*)/* (z')f^(z'). 



Projecting this onto the bubble wall zero mode then transforms Khlebnikov's result (ET 
into: 



a, b 



U 



dzfr(z)d z M 2 (Mz))f*(z) 



• (2-12) 



This will reproduce the result from (|1.8|) provided the sum over (a, b) above can be appro- 
priately related to the reflection and transmission coefficients 1Z and T. 



Using the eigenequation (|2.5| ) that defines the it is straightforward to show that 
the ^-integrand is a total derivative: 



f%*(z)d z M\<t> cl (z))f b R (z) = d z 



1 <->2 

-f a *d f b 



(2-13) 



where d z is defined by 



fd z g = (d 2 J)g - 2(d z f)(d z g) + f(d 2 z g). 

*2 



(2-14) 



So the result ( |2. 12| ) depends only on the values of /"* d z y K at infinity, where becomes a 
superposition of incoming and outgoing plane waves. On plane waves, note that d z gives 



(e lkz Y d z e lkz = -4k 2 , {e~ lkz Y 8 z e lkz = 0. 



(2-15) 



Now let us choose the basis for fj* depicted in fig. [5] and parameterize the asymptotic 
behavior as 

^ (e^ z + S ++ e~^ z ) , s--oo, 

(2-16) 



OO , 



where the 5 are functions of k and 



i ( e - ifc »* + 5__e ifc3 *) , x^+oo 



k 2 = k\ + M 2 (-oo) = fc§ + M 2 (+oo) 



(2.17) 



The matrix of coefficients 



5++ S+. 



(2.18) 



is symmetric and unitary.* The reflection and transmission coefficients are 1Z = l^-i-l 2 = 

IS 1 --! 2 and T = | | 2 = \S y\ 2 . The normalization condition (|2.6|) is verified in the 

appendix. (When n 2 is too small to allow transmission from phase 1 to phase 2, then 
1Z = 1, T = 0, and the /+ solution should be dropped.) 



Using the eigenfunctions (|2.16|) , the properties of 5 1 , and the action ( |2.15|) of d z , it is 
straightforward to find that 



E 

a, b 



1 ^2 , 

_ fa* o rb\+°° 

2 J ft " z J k | — oo 



1 



2tt 2 



[2{k 2 + kj)TZ + (h - k 2 ) 2 T] 



(2.19) 



[When k 2 is too small to allow transmission from phase 1 to 2, this instead becomes simply 
(2/c 2 )/2-7r 2 .] Now use this, together with (|2.13| ), to rewrite the result ( |2.12j ) for the bubble 
wall velocity in terms of 1Z and T. As promised, the result is identical to the result (|1.8| ) 
of the last section, where I reviewed the simple picture of friction caused by reflection and 
transmission at the bubble wall. 



3. A Toy Example 

I now want to specialize to Khlebnikov's toy example and understand his results in 
terms of reflection and transmission probabilities. Khlebnikovit took M 2 (<fi) = Mq + 2gcj) 
in the limit AM 2 < Ml < T 2 where AM 2 = M 2 (z= + oo) - M 2 (z= - oo). Consider first 
the transmission term in the net wall pressure ( |1.8|) , and consider the contribution to the 
integral from momenta AM 2 <C p 2 -C Mq. T is approximately 1 for such momenta, and 
it is easy to extract the behavior of the integrand in these limits. Taking E 2 ~ Mq 2 + p^_, 
Piz ~ P2z, Piz — P2z ~ AM 2 /2p z , and n(E) ~ T/E, this contribution to the transmission 

* This follows from examining the Wronskians W of / + with f~ or f~* and equating 
W(z= - oo) with W(z= + oo). 
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term in (|1.8| ) becomes 

12 



0U I ^0f n ^ n ^ + 1 ^-P^ 2T6 ^ 



u(AM 2 ) 2 T-^-r / ^ / a ^ (3.1) 



4(27r) 3 i p 2 7 (^+M 2 )3/2 
(AM 2 ) 2 T / Mq \ 
W W 327r 2 M V AM 2 J ' 

where p z has been integrated over the region AM 2 <C p 2 <C Mq, and where I have ignored 

terms on the right-hand side that are not enhanced by a logarithm. It is easy to check 

that this contribution dominates (by the logarithm) over the other regions of integration, 

p 2 <AM 2 and M 2 <p 2 . 

Now consider the reflection term in the net wall pressure ( |1.8| ). Generically, 1Z falls 

to zero for p 2 3> AM 2 , and so one might not expect to get a logarithmic enhancement as 

in ( |3.1| ). However, if the bubble wall were a step function, with zero width, then 1Z falls 

slowly enough at large p z that one still gets a logarithm. Specifically, for a step potential, 

Piz ~ P2z 2 



K= Flz F2z . (3.2) 

PU + P2 Z 

Considering the region AM 2 <C p 2 <C Mq as before (so that p lz ~ p 2z and T ~ 1), then 

2(p 2 z +p 2 2z )n^(p lz -p 2z ) 2 T, (3.3) 

and one sees that the 7Z and T terms in (|1.8|) contribute equally, each giving (|3.1|). If 
the spatial width 5 of the wall is small but non-zero, then 1Z will behave like ( |3.2|) for 
p z <C 1/5 but will fall faster with p z than ( |3.2| ) for p z ^> 1/5. The logarithm found for the 
transmission term in ( |3.1| ) will then, for the reflection term, be cut off by 1/5 rather than 
Mq if 1/5 < Mq. At the level of the leading logarithm, the various cases are 



I3u J ^J^En(E)[n(E) + l]2(p? z + P L)RS(fi, 



W,,i ^7r- < It, (£, 



\ ln U^) ' am 2 <5- 2 <m 2 , 

0, <T 2 <AM 2 , 
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where zero in the last case means that there is no term with a logarithmic enhancement. 
Putting the approximations ( |3.1|) and ( |3.4j) into the formula ( |1.8|) for the net pressure, the 
result for the wall velocity in Khlebnikov's model is 

U= {AM 2 ) 2 Tl + ° {U/1) > (3 - 5) 



where 



I 



4^)+l*GSr)> AM*<5-*<Ml (3.6) 
In (A) . 5~ 2 <AM 2 . 



For a specific form of c i(z), one could also calculate the non- logarithmic correction to / 
by precisely computing the reflection and transmission coefficients. 

Khlebnikov gave a more approximate analysis of this model where he (1) accidentally 
missed the transmission term,* and (2) computed the contribution of fig. [|b to the ITr 
in the approximation that the background field (p is ignored, so that the $ propagator is 
simply l/(-P 2 — Mq ). This approximation turns out to be equivalent to taking AM 2 — > in 
the logarithms ( |3.6j ) and produces a logarithmic infrared divergence. He was forced to cut 
off this divergence by including rescattering effects, which come from resumming higher 
order of the loop expansion. However, we see now that there is indeed no logarithmic 
divergence in the simple one- loop result. 

In the case of the electroweak phase transition, Mq and AM 2 are not separate scales, 
and so no logarithmic enhancement of the friction should be expected. In electroweak 
models, the wall width is typically large compared to inverse particle masses, and then the 



* The intended step between eqs. (19) and (20) of ref. || requires the identity 

S(E - Ex) = (cos9- + -j^rAk), 
pk \ 2p J p\ cos 6*1 

but the last term was left out. There is also a factor of 2 mistake going from eq. (20) to (21) 
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reflection and transmission coefficients can be well approximated by simple theta functions, 
0{—P2 Z ) and 0(p2 Z ), representing complete transmission when energetically allowed. This 
approximation was used in refs. 

I have shown that the one-loop application of Khlebnikov's more general fluctuation- 
dissipation result is equivalent to the simple physical picture of computing the friction on 
the wall from particle collisions in the ideal gas limit. For the electroweak phase transition, 
however, one generally may not ignore the effects of rescattering, and the reader should 
consult refs. fl|-|3[] for attempts to analyze rescattering effects. 

I thank Lowell Brown for many long and useful discussions, for reminding me how to 
do quantum mechanics, and for getting me interested in this problem. I also thank Larry 
Yaffe and Sergei Khlebnikov for useful discussions. 



Appendix A. 



In this appendix, I check that the solutions (|2.16| ) are normalized according to (| 
For a complete set of states, this normalization condition is equivalent to 



dK 2 fr(z)f:(z f ) =5(z- z 1 ). 



(A.1) 



The normalization of the basis (|2.16|) may be checked by verifying (|A.1[) for any particular 
range of z. Let's focus then on z, z' — > — oo. Using ( |2.16| ) and the properties of S, one 
finds 



i 



Aixk A 



ik 1 (z-z') _|_ e ik 1 (z-z') _|_ 2S++e~ ikl ( 2-1 " 2 ') -f 2S^_ Jr e ikl ( zJrZ '^ 



(A.2) 



Then, since k 2 — > fcf , 

J d K 2 r K *(z)f a K (z') - S(z ~z') + l dh[S ++ + S* ++ ]e- ik ^ z+ *'\ (A.3) 
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The last term vanishes due to two analytic properties of S ++ : (1) S* r+ (ki) = S ++ (—ki), 
and (2) S++(k\) is analytic in the upper half-plane provided the potential has no bound 
states (which, in applications of interest, is generally the case). The last term can therefore 
be written as an integral from — oo to +00, which vanishes by closing the contour in the 
upper-half plane. The first property follows from the fact that, if fi^(z) is a solution, 
then f£ (z) is also a solution. For insight into the second property, suppose that S++(ki) 
has a pole in the upper-half plane. At this pole, k^ = \/k\ — AM 2 will also be in the 
upper- half plane. Now renormalize the solution (|2.16|) for /+ by dividing it by <S++, and 
then take ki to be the position of the pole so that S++ — > 00. The "incident" wave term 
disappears. The solution falls exponentially as x — > ±00 and therefore describes a bound- 
state, which contradicts the assumption that there are none. (If there are bound states, 
the contribution from the poles of S++ reflects the fact that the sum over bound states 
should have been included on the left side of ( |A.1| ).) The simple argument above doesn't 
prove that there aren't other types of singularities in the upper-half plane, and the reader 
should try ref. [[nj for a more general discussion of the analytic properties of the S matrix. 
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Figure Captions 

Fig. 1. Qualitative form of the effective, finite-temperature Higgs potential during a first- 
order electroweak phase transition. 

Fig. 2. Qualitative profile of the bubble wall. 

Fig. 3. Schematic depiction of two solutions to the Schrodinger equation describing a 
particle propagating in the background of the wall. A wave may be incident from 
either (a) the left or (b) the right, and it reflects and transmits at the wall. The 
heavy solid line depicts the potential. 

Fig. 4. One-loop contributions of $ to the Higgs self-energy. 
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